Analogies Between Scaling in Turbulence, Field Theory and Critical Phenomena 
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We discuss two distinct analogies between turbulence and field theory. In one analogue, the 
field theory has an infrared attractive renormalization-group fixed point and corresponds to critical 
phenomena. In the other analogue, the field theory has an ultraviolet attractive fixed point, as in 
quantum chromodynamics. 
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^^ ■ In this brief note, we shall discuss two distinct analogies between turbulence at high Reynolds number, statistical 

Q\ \ field theory and critical phenomena. Such analogies are implicit in attempts to make a useful connection between 

these apparently disparate phenomena, but a number of points deserve to be made explicitly. Some of our remarks 

are essentially trivial observations; nevertheless their ramifications have not always been respected by approximate 

theories §J. 

The fully-developed turbulent regime is specified by the condition that the integral length-scale L of the largest 
eddies should be much larger than the length scale rjd at which dissipation is effective. The latter is usually defined 
as the inverse of the wavenumber kd at which occurs the peak of the enstrophy spectrum k 2 E{k). Equivalently, the 
Reynolds number Re = U L/v should be large, where U is a characteristic large-scale velocity and v is the molecular 
viscosity. We consider the statistical steady-state with constant mean rate e of energy injection by the turbulence 
production mechanism, which, by stationarity, is also the mean rate of energy dissipation by molecular viscosity. 
We shall take the point of view that the problem of turbulence is to understand the properties of the stationary 
(but presumably non-Gibbsian) probability distribution governing the velocity field of a turbulent fluid; it is this 
distribution which is used below when we write averages in the context of turbulence. 

The first analogue of turbulence is to field theories with an infra-red attractive renormalization-group (RG) fixed 
point. This is the situation for field theories corresponding to critical lattice spin-systems. For example, consider 
C^ ■ a system of spins <r(r) in d spatial dimensions, where r runs over a lattice aZ d of lattice-constant a, governed by a 
(dimcnsionlcss) Hamiltonian of the form: 
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with e^, fi = 1, 2, . . . , d being unit lattice vectors. The parameter r depends upon the temperature of the spin system, 
with an order-disorder transition occuring at the critical value r c . In the critical regime (r — t c )/t c <C 1 there is, for 
large space-separation, universality from the specific short- distance, or lattice-scale interactions. On the other hand, 
turbulence scaling in the high wavenumber inertial-range is believed to exhibit universality from the small wavenumber 
stirring mechanisms. Therefore, in this analogy the roles of space and wavenumber are interchanged. A detailed list 
of correspondences (and definitions) may be drawn up for an analogy along these lines, which is motivated below: 

Turbulence Critical Phenomena 

space separation r wavenumber k 

viscosity, v temperature variable t — t c 

energetic length-scale, L UV cut-off A, (or, inverse lattice-spacing, a -1 ) 

mean dissipation, e stiffness constant, K 

dissipation wavenumber, kd = r]J correlation length £ 

velocity correlation function. spin correlation function, 

5 2 (r) = <(v(r' +r) - v(r')) 2 ) C(k) =E reaZi a d e^(a(r)a(0)) 

intcrmittency exponent, /i correlation exponent, rj c . 

This analogue is similar to that pointed out previously by Nelkin g, deGennes g and Rose and Sulem J4|, the 
critical (r — ► r c ) limit of equilibrium spin systems being considered analogous to the zero- viscosity limit of turbulence. 
However, we have made a correspondence of the constants e and K which does not seem to have been pointed out 



before. Also, we propose that \i and r\ c are directly analogous, rather than being related as fj, <-> | — r) c: as suggested 
by the earlier authors. 

The above analogy can be motivated by comparing Kolmogorov's 1941 theory of turbulence |j| and Landau's 1937 
"mean-field" theory of critical phenomena J6[. In both theories, critical exponents are obtained by dimensional analysis 

In the Landau mean-field analysis, it is implicitly assumed that the long wavelength properties are independent of the 
lattice constant a. Thus the limit a — + of spin-correlations is presumed to exist, an assumption that the asymptotics 
are of the first kind J8|. In that case the only remaining parameter with units of length is r = ((r — r c ) / ' K) , 

— 1/2 

yielding the Landau prediction for the correlation length £l ~ r . The values for the critical exponents follow from 
this assumption: for example, 

C(k) ~K-k~ 2 , (2) 

for kt; 3> 1, where it has also been assumed that the critical limit ro — * exists. In fact, Landau's first assumption 
that the limit a — ► exists is generally wrong, and the correct asymptotic scaling behavior is of the second kind [B 

C(k)~K ■k- 2 {ka) ri ", (3) 

for some r\ c > 0, in which the microscopic length a appears explicitly. Likewise, the Landau scaling of the correlation 
length has an "anomalous" correction arising from a-dependence: £ ~ £l • (roa 2 ) for some 9 > 0, or £ ~ Tq V ° with 

Vc = \-e. 

In the same way, Kolmogorov assumed that the limit L — > +oo should exist in turbulence with finite velocity 
correlations. In that case, the only remaining length-scale is the dissipative scale estimated as r\K ~ {e)~ 1 / A u 3 / i . 
Thus, Kolmogorov obtained by dimensional analysis 

S 2 (r) ~ (er) 2/3 , (4) 

when kdT 3> 1, if the limit v — ► exists for velocity correlations. As with the corresponding assumption in Landau's 
theory, Kolmogorov's hypothesis of existence of the first L — » +oo limit is questionable, due to the build-up of 
fluctuations of energy flux in the energy cascade. (Ironically, this criticism originates in part from a remark of Landau 
himself |l!). Instead, simple cascade models indicate that the scaling law should instead be of the form 

ft(r) ~ (er)^ (?) \ (5) 

for some /i > 0. Similarly, these models suggest p0[ that the dissipation scale may have a slight dependence upon L, 
as r)d ~ r\K ( J*- ) for 5 ^ 0, leading to a scaling r]d ~ rfc with w = 1 — 5. 

It is clear that in these formulas, K,a,r] c ,v c are homologous, respectively, to e, L -1 ,/i,u;. Notice that typically ry c 
is small in critical systems even in two or three dimensions and that /i represents an empirically small modification 
to the dimensional analysis result of 5/3 for the exponent appearing in the Fourier transform of 52: 

E(k) - (e) 2/3 k- 5/3 {kL)-». (6) 

The constants K and e play a similar formal role in the two theories. To see this, consider a Martin-Siggia-Rose 
(MSR) Lagrangian fll| , |l^| for the steady-state turbulent cascade state produced by driving with a Gaussian random 
force f , white noise in time, with zero mean and covariance 

(/i(r,*)/ i (r' ) 0> = 2%F(r - r')5(t f). (7) 

It is easy to check that F has the units of energy dissipation and, indeed, it is directly related to the mean energy 
injection rate by e = 2F(0). Therefore, the quadratic term in the MSR action corresponding to the "response field" 
v is 

S(°)[v] = i I dt f drdr'v(r,t)F(r-r')v(r',t), (8) 

which is proportional to e. S^ has a formally similar structure to the quadratic term in the spin-Hamiltonian, which 
is proportional to K. 



A second, rather different analogy can be made between turbulence and field theory with an ultraviolet (UV) 
attractive fixed point, or, field theory for short. (It is generally accepted that, in order for a continuum field theory 
model to be well-defined, it must be "asymptotically safe," i.e. it must correspond to a lattice model with a UV 
attracting fixed point.) Field theory, such as UV asymptotically free quantum chromodynamics (QCD), exhibits 
scaling at short distances, just as turbulence is believed to do. Therefore, in this analogy space corresponds to space. 
The spin system we examined before may still be used in this analogy, when it is considered in dimension d < 4 along 
the single RG trajectory which flows in the UV direction into the non-Gaussian Wilson-Fisher fixed point (previously 
we considered the theory at a generic point on or slightly away from the critical surface). It is convenient to use the 
field-theoretic notation, introducing a (lattice) field 

0(r) = K^Mr), (9) 

which has dimension ^(d — 2) in units of inverse-length. As before, a detailed list of correspondences may be drawn 
up: 

Turbulence Field Theory 

space separation r space separation r 
viscosity v (or, Kolmogorov scale, t\k) lattice-spacing a 

energetic length-scale, L correlation length £ 

Kolmogorov wavenumbcr, k n UV cut-off A = a -1 

velocity correlation function, field-theoretic Green function, 

S 2 (r) = ((v(r' + r) - v(r')) 2 ) G(r) = (0(r)0(O)> 

These correspondences essentially all follow from the comparison in the last line. For example, the length-scale L gets 
its name "integral-scale" from the fact that it is an integral correlation length of the longitudinal velocity correlation 

f(r) = ((f • [v(V + r) - v(r')]) 2 ), defined by 

Therefore, it corresponds naturally to the correlation length £ defined by the exponential decay rate of G(r). Likewise, 
the dissipation length tjk provides a short-distance cut-off for turbulent velocity correlations, analogous to the lattice 
cutoff a for the field theory correlations. The correspondences made here, (k v ,L) «-► (A,£), are the opposite of 
those made in the first analogy, (k Vl L) <-► (£, A). In the present analogy, the zero- viscosity limit is analogous to the 
continuum limit a — > of field theory. 

Having pointed out the analogies between turbulence and field theory, it is appropriate to emphasize certain 
important distinctions. First, we discuss the issue of universality. In turbulence, short-distance scaling is believed to 
be generic and essentially independent of large-scale statistics or driving mechanisms. Scaling laws, strictly speaking, 
require an inertial range of infinite extent, so that we confine our discussion here to the idealized situation of zero 
molecular viscosity (i.e. the critical limit according to our first analogy). We will assume that the stationary probability 
distribution governing the turbulent state is a fixed point in a function space, whose axes may be thought of as coupling 
constants for all conceivable local (and perhaps non-local) operators. This distribution is usually assumed to be the 
fixed point of a RG transformation which integrates out short wavelength degrees of freedom, just as in critical 
phenomena. (An alternative to this procedure will be mentioned below.) Short-distance universality implies that the 
ultra-violet RG flows should all be in towards the fixed point: this distribution is a global sink in the ultra-violet. (Here 
again we ignore crossover phenomena associated with finite molecular viscosity or other short-distance regularization 
mechanisms.) Thus, the usual infra-red RG flow diagram, moving the system to longer and longer length scales, will 
contain a fixed point with a presumably infinite number of relevant directions. The physical significance of relevant 
directions, in the context of critical phenomena, is that they represent the parameters that must be tuned in order for 
the system to be at criticality (e.g. temperature is at the critical temperature, external field is zero). In the context 
of turbulence, these relevant directions correspond to the myriad of different large-scale stirring mechanisms which 
generate the same short-distant behaviour. One of these relevant directions, corresponding to the temperature, is the 
"eddy viscosity" , which is the effective viscosity at a given scale generated by turbulent degrees of freedom at shorter 
scales. The eddy viscosity will tend to the bare, molecular viscosity (assumed zero) as the turbulence fixed-point is 
approached at short length-scales. In contrast, even field-theories with a UV fixed point must generally be tuned to 
lie on the low-dimensional surface which attracts to the fixed point, and as mentioned above, the number of variables 
that must be tuned is equal to the number of IR unstable directions at the fixed point: this is usually a finite number. 
This general argument does not distinguish between relevant and marginally relevant directions. 



The second distinction is that scaling behavior may be quite different in turbulence and typical field-theories. The 
latter usually show simple "gap-scaling " or "hyperscaling" in which higher-order correlations scale with exponents 
which are just integer multiples of a single field dimension. For example, if d$ is the scaling dimension of the field <fi, 
then typically 

(<K*ri)---<KXr P ))~\- p - d *, p= 1,2,3,... (11) 

as A — > 0. In contrast, higher-order velocity structure functions, S p (r), are believed to scale in high Reynolds number 
turbulence with r — > as 

<(f • [v(r' + r) - v(r')]) p > ~ {erf'* (?) " , (12) 

for x p < 1 13 13]. This type of "multifractal" scaling is found in very few field theories. In the cases where it occurs, it 
is usually associated with sequences of field variables with negative scaling dimensions. More generally, the multifractal 
scaling occurs if there is a sequence of variables "additively coupled" in the operator product expansion and negative 
scaling dimensions are one common mechanism by which such sequences are produced |lq,Eq] . A representative case is 
the 0(N) nonlinear er-model in dimension d — 2+e with N < 2 flI7f| , which has such a sequence. This example is closely 
related to Wegner's theory of the mobility edge in disordered electron systems and the hypothesis there of multifractal 



electron wavefunctions for the localized states 17 1. Note that any variable with scaling dimension x < d corresponds 
to an unstable direction at the fixed point. Therefore, the variables with negative scaling dimensions correspond to 
an infinite set of unstable directions, and the associated fixed point is UV attractive in a large domain, as suggested 
in the preceding paragraph. Motivated by this type of example, an RG theory of high Reynolds number scaling in 
turbulence can be developed [ |19| , p0[ in which the "anomalous dimensions" are the negative scaling dimensions x p 
of the powers of the velocity-gradients, (Vv) p . In traditional turbulence terminology this corresponds to generalized 
flatnesses for velocity gradients diverging as powers of the Reynolds number, which is also observed in both simulation 
and experiment flq , pifl . The same phenomena lead also to the possibility of a hierarchy of dissipation lengths rff 
delimiting the short-distance end of the scaling range of the pth structure functions S pi each having a different 

dependence on molecular viscosity ry£' ~ r]^ F) [jl(| . 

The third distinction that we wish to mention is the form of the continuum limit. In field theory taking the 
continuum limit requires making a multiplicative renormalization of the lattice field- variables, 

0](r) = Z{a) ■ 4>{t), (13) 

where Z(a) ~ a~ 7 * and 7^ is an "anomalous dimension" of the renormalized field [</>]. The necessity of this renormal- 
ization is connected to the UV divergences in the field-theoretic Greens functions of the "bare" fields (cf. Eq.(ph). 
Unlike the fields [<fi] which exist only as distributions in the limit a — > 0, the velocity field in turbulence must remain an 
ordinary function in the limit v — > 0. The requirement of finite mean kinetic energy {\v 2 ) < 00 implies this (since the 
velocity field must then be locally square- integrable with probability one). On the other hand, the velocity gradients 
will exist necessarily only as distributions and UV divergences may appear associated to their products at a single 
space-point. Only these products require "renormalization" (by suitable powers of t\k ) in the zero- viscosity limit 
and any "anomalous dimensions" in inertial-range turbulence scaling must therefore be associated with the velocity 
gradients rather than the velocities themselves. 

The two analogies we have pointed out differ most obviously in the reverse roles played by scale (wavenumber) and 
space. However, they differ more essentially in terms of their motivation, the first (with critical phenomena) being 
more suggestive in physical terms and the second (with field theory) arising from mathematical considerations of the 
formal renormalization procedure. 

In our first analogy, the correlation length in critical phenomena and the dissipation scale in turbulence are intrinsic 
properties of the system, determined by its detailed dynamical and statistical properties, rather than fixed external 
inputs. In contrast, the lattice spacing in the first case and the integral scale in the second case are parameters fixed 
by the experimental setup or the definition of the model. The inverse role of large and small scales in the two cases 
arises from the different character of "cascade" in the two instances. According to the cascade picture, there is a 
transfer of excitation on the average from the large turbulent eddies to the small ones by a stepwise process, which 
is chaotic in nature and entails a loss of memory of the large-scale statistics. Wilson has emphasized p2fl that there 
is also a "cascade of fluctuations" in critical phenomena. Droplet fluctuations nucleated at the lattice-scale in the 
critical state can grow to the size of the correlation length (and vice versa). However, it is now the details of the 
lattice structure which are lost and the scale-invariant distributions of the large "droplets" which are universal. 

These facts have suggested to several authors p3| , p3j , |25J that in constructing an RG theory of turbulence, it may 
be more natural to reverse the usual procedure and to eliminate low-wavenumber shells rather than high- wavenumber 



shells as in the Wilson method. Some cautionary remarks on such a procedure at a physical level are made in p6[ . 
However, our discussion of the field theoretic analogy shows also that to study short-distance UV-scaling, an IR- 
climination RG is not required. Lattice QCD is a good case in point, where the short- distance asymptotic freedom 
is studied by the same UV-elimination procedure as used in critical phenomena. The UV-elimination RG is also 
closely related to an important practical problem of turbulence theory, that of "subgrid-scale eddy modelling," and is 
naturally described in terms of simple "eddy viscosity" concepts. Indeed, the UV-elimination RG is deeply connected 
with such eddy viscosity ideas, and RG-invariance under this operation is nothing more than a restatement of the 
fact that the inertial range theory is independent of the molecular viscosity [£(J . 
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